FTUV/96-29 
IFIC/96-33 
HEPHY-PUB 647/96 
UWThPh-1996-36 

New signatures for a light stop at LEP2 in SUSY 
models with spontaneously broken R-parity 

A. BartftW. Porodi 

Institut fur Theoretische Physik, University of Vienna 
A- 1090 Vienna, Austria 

3 ■ 

M. A. Garcla-Jareno I, M. B. Magroi and J. W. F. Valle I 

Instituto de Fisica Corpuscular - C.S.I.C. 
Departament de Fisica Tebrica, Universitat de Valencia 
in ■ 46100 Burjassot, Valencia, SPAIN 

\q ; yittp: 1 1 neutrinos.uv.es 



o. 

)0- W. MajerottoH 



Institut fur Hochenergiephysik, Akademie der Wissenschaften 
A- 1050 Vienna, Austria 



Abstract 



X 

In a class of supersymmetric models with R-parity breaking the 
lightest stop can have new decay modes into third generation fermions, 
t\ — > 6 + t. We show that this decay may be dominant or at least com- 
parable to the ordinary R-parity conserving mode ti — >• c + x±, where 
Xi denotes the lightest neutralino. The new R-parity violating decay 
modes could provide new signatures for stop production at LEP. 
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1 Introduction 



So far most searches for supersymmetric particles have so far been performed 
within the framework of the Minimal Supersymmetric Standard Model 
(MSSM) [JTp. Although the MSSM is by far the most well studied realization 
of supersymmetry, there is considerable theoretical and phenomenological in- 
terest in studying possible implications of alternative scenarios [g], in which 
R-parity is broken. Indeed, neither gauge invariance nor supersymmetry re- 
quires the conservation of R-parity. The violation of R-parity could arise explic- 
itly [Q] as a residual effect of some larger unified theory [g], or spontaneously, 
through nonzero vacuum expectation values (VEV's) for scalar neutrinos []5p. In 
realistic spontaneous R-parity breaking models [§, [7|, |3p the scale of R-parity 
violation lies around the TeV scale. Its effects can therefore be large enough 
to be experimentally observable. Moreover, these models are fully consistent 
with astrophysics and cosmology [§]. In either case the supersymmetric (SUSY) 
particles need not be produced only in pairs, and the lightest of them could decay. 

There are two generic cases of spontaneous R-parity breaking models to 
consider. If lepton number is part of the gauge symmetry there is an additional 
gauge boson which gets mass via the Higgs mechanism, and there is no physical 
Goldstone boson In this model the lightest SUSY particle (LSP) is in general 
a neutralino which decays mostly into visible states, therefore breaking R-parity. 
The main decay modes are three-body decays such as 

Xi - ffv, (1) 
where f denotes a charged fermion. Its invisible decay modes are in the channel 

Xi - 3za (2) 

Alternatively, if spontaneous R-parity violation occurs in the absence of any ad- 
ditional gauge symmetry, it leads to the existence of a physical massless Nambu- 
Goldstone boson, called majoron (J). Thus in this case the lightest SUSY parti- 
cle is the majoron which is massless and therefore stable 0. As a consequence 
the lightest neutralino Xi m ay decay invisibly as 

xl^v + J, (3) 

7 The majoron may have a small mass due to explicit breaking effects at the Planck scale. 
In this case it may decay into neutrinos and photons. However, the time scales are only of 
cosmological interest and do not change the signal expected in laboratory experiments JlOj]. 
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where the majoron is mainly a singlet [7p. This last decay conserves R-parity 
since the majoron has a large R-odd singlet sneutrino component. 

In this paper we focus on the decay modes of the lightest top squark in 
models where supersymmetry is realized with spontaneous R-parity violation. In 
such models the lightest stop could even be the lightest supersymmetric particle 
and be produced at LER Indeed present e + e~ collider data QUI ] as well as pp 
collider data QI^ ] do not preclude this possibility. 

In supersymmetric models with spontaneous breaking of R-parity or in 
models, where this violation is parametrized explicitly through a bilinear super- 
potential term of the type £H U Ql3|], the stop can have new decay modes such 
as 



due to mixing between charged leptons and charginos. In this paper we show 
that this decay may be dominant or at least comparable to the ordinary R-parity 
conserving mode 



where Xi denotes the lightest neutralino. 

2 Lepton-Gaugino-Higgsino Mixing 

The basic tools in our subsequent discussion are the chargino and neutralino 
mass matrices. The chargino mass matrix may be written as [|6p 



b + r 



(4) 




(5) 



H+ -iW+ 




heijVd -KijV Rj V2g 2 v Li 
-h e ijVu /i V2g 2 v d 
V2g 2 v u M 2 



(6) 



Its diagonalization requires two unitary matrices U and V 



(7) 



(8) 



where the indices i and j run from 1 to 5 and ipf = (e|, e\ , e% , H+, — iW + ) and 
^7 = ( e i > e 2 > e 3 ' ^7> -iW~). 
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While the details of the neutralino mass matrix are rather model dependent, 
for our purposes it will be sufficient to use the following effective form given by 







H u 


H d 


-iW 3 


-iB 












givu 


-givu 


H u 


hi/ijVjlj 





-fi 


-giv u 


giv u 


H d 





-V 





92V d 


-giVd 


-iW 3 






giVd 


M 2 





-iB 


-givu 


giv u 


-giv d 





Mi 



(9) 



This matrix is diagonalised by a 7 x 7 unitary matrix N, 

X° = (10) 

where ip® = (ui,H u ,H d , — iW 3 , —iB), with z/j denoting the three weak-eigenstate 
neutrinos. 

In the above equations v u and v d are the VEV's responsible for the breaking 
of the electroweak symmetry and the generation of fermion masses, where the 
combination v 2 = v\ + v\ is fixed by the W, Z masses. v R is the VEV mostly 
responsible for the spontaneous violation of R-parity i. Moreover, Mi 2 denote the 
supersymmetry breaking gaugino mass parameters and g^ 2 are the SU (2) g) U{1) 
gauge couplings divided by \/2, and we assume the canonical relation M 1 /M 2 = 
^tan 2 8 w . Note that the effective Higgsino mixing parameter fi may be given in 
some models as ^ — ho (}, where () is the VEV of an appropriate singlet scalar. 

Most of our subsequent analysis will be general enough to cover a wide 
class of SU(2) ® U(l) models with spontaneously broken R-parity, such as those 
of ref. [||, Q7p, as well as models where the majoron is absent due to an enlarged 
gauge structure []9p. Many of the phenomenological features relevant for the LEP 
studies discussed here already emerge in an effective model, where the violation 
of R-parity is introduced explicitly through a bilinear superpotential term of the 
type IH U [ITSp. In this case the combination h v v R is replaced by a mass parameter 



€. 



Typical values of interest for LEP2 for the SUSY parameters /i, M 2 , and the 



8 There is also a small seed of R-parity breaking in the doublet sector, vl = {vlt), whose 
magnitude is related to the Yukawa coupling h u . Since this vanishes as h v — > 0, we can naturally 
obey the limits from stellar energy loss 
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(11) 



parameters h vi $ lie in the range given by 

-1000 GeV < n < 1000 GeV 40 GeV < M 2 < 200 GeV 

10 10 < h„ 13 , h v23 < 10 1 10- 5 < h u33 < 10 1 

while the expectation values can be chosen as 

vl = v L3 = 100 MeV v L1 = v L2 = 

50 GeV < v R = v R3 < 1000 GeV v R1 = v R2 = 

with 1 < tan/3 = ^ < 40. 

There are restrictions on these parameters that follow from searches for 
SUSY particles at LEP Q and the TEVATRON In addition, we take 



(12) 



into account the constraints from neutrino physics and weak interactions phe- 
nomenology []15[j, which are more characteristic of R-parity breaking models. 
These are important, as they exclude many parameter choices that are other- 
wise allowed by the constraints from the collider data, while the converse is not 
true. Due to these constraints R-parity violation effects manifest themselves 
mainly in the third generation as expressed in eq. (|I2"|). In order to evaluate 



the masses and couplings of charginos and neutralinos we have performed a 
sampling of parameters in our model which are allowed by all constraints above. 

One finds that the R-parity conserving couplings for the lightest neutralino 
Xi and the lightest chargino xt are of the same order as those in the MSSM. The 
R-parity violating chargino couplings may reach a few percent or so for chargino 
masses accessible at LEP [ ]T5j , |IEp. These couplings will induce a mixing between 
the charginos and the tau lepton, leading to the R-parity violating decay in eq. (g). 
The corresponding decay width will be calculated below. 



3 Decays of the Top Squark 

As we saw in the previous section, R-parity breaking effects should occur mainly 
in the third generation, due to the mixing of the charged leptons with the 
charginos. As a result the decay mode ti — » b + r is possible. In the following 
we assume < m b + m -+ . In this case the most important R-parity conserving 
stop decay at LEP2 is t\ — >• c + Xi Now the question arises whether the 

branching ratios of the decays into br and cx\ are of comparable size. As we will 
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show ti — > b + r can be equally or even more important than t\ — >■ c + x°- We 
have also taken into account the decay t\ — > c + X2, which is also possible in a 
small region of parameter space. In the alternative case mi > nib + the 
decay mode t± — > 6 + Xi~ would clearly dominate the decays (f|) and @. A detailed 
discussion of all R-parity conserving stop decays at LEP2 is given in yi£j[J. 

The decay ?i — > c + x? is only possible in higher order (see Fig. p. As 
shown in [|15p the dominant contribution is due to ti — > c L — > cx?. This can 
be parametrized by the following squark mass squared matrix in the (t L , t R , c L ) 
basis: 



Explicit expressions for Mf L , Mf ,a t m t , M? L , L and R can be found in ref. [jlffi. 
Quite generally, flavour changing mass terms are induced in supergravity mod- 
els assuming unification and evolution according to the renormalization group 
equations Q20|]. The light stop eigenstate can therefore be written as: 



Estimates for 5 were given in p^[]. These are, however, very model dependent as 
pointed out in U21|j. Nevertheless a value 5 ~ y^t m ay be regarded as typical. 
For our present phenomenological discussion we take a very strong t-c mixing, 
S ~ 0(0. 1). To be specific, we shall take 5 = 0.1. Note that this is a conservative 
assumption, i.e. should 5 be smaller, the relative importance of our novel decay 
eq. @ would be further enhanced. 

The decay widths for (£|) and @ are: 




(13) 



h = c x t L + c 2 t R + 5c L . 



(14) 



(fi -> b + r) 




x((/ 2 + k 2 )(m 2 {i - mf - m 2 .) - 4lkm b m T ) 



(15) 



(ii^c + x°) 




la 1 



(16) 



where 





U u ci (17) 



/i 



(tan 9 W N 47 + 3N 46 ), f 2 = - ^— (tan ^iV 57 + 3iV 56 ) 



(18) 
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Here V34, U34 (V35) denote the amount of tau-Higgsino (tau-gaugino) mixing (see 
eq. @ and eq. @) and K tb = 0.999, K cb = 0.04 are the corresponding CKM matrix 
elements [1221]. 



For a dominantly Higgsino-like neutralino []23|] (small N i6 and N i7 ) the R- 
parity conserving decay is induced by the charm Yukawa coupling and is there- 
fore suppressed (see eq. ( |16D and eq. (|18|)). As a result, the branching ratio of the 
R-parity violating br mode will be enhanced. 

In our present numerical analysis we have fixed m tl = 80 GeV, c-i = 0.5, v R = 
100 GeV, vl = 0.1 GeV and ^33 = 0.03. We have varied the other parameters 
in the range given in eq. ( |12| ) and eq. ([11]). As already mentioned, we have 



taken into account all bounds from the LEP experiments as well as neutrino 
physics []T^] I. In Fig. 2 we show the maximum value of the branching ratio for 
ti — > b + t in per cent as a function of the lightest neutralino mass m^o in the 
range 1 < tan (3 < 40. One can see that BR(fx — > b + r) can easily reach 80 % for 
neutralino masses of 50 GeV. There are a few points, not shown in Fig. 2, which 
lie above the "maximum curve" even for relatively small neutralino masses. 
They correspond to "fine-tuned" situations where f± almost vanishes due to a 
cancellation between the terms in eq. (|T8|) for small tan /3 values. 

Particularly interesting is the dependence on tan (3, which can be seen in 
Fig. |3[ The most characteristic feature is the existence of a strong correlation 
between the branching ratio for the br stop decay mode and tan/5, apart from a 
small region in parameter space characterized by small tan j3 where fi vanishes 
and the R-parity conserving decay is suppressed 0. One can clearly see that 
most points with a large branching ratio of £1 — > b + r have a large tan j3. Note 
that in the framework of the MSSM with radiative electroweak breaking such 
large tan/3 models have been advocated |]25[j. The trend towards large tan/3 can 



be understood as being due to the influence of the large bottom Yukawa coupling 
(see eq. (|I7|)). It is also interesting to note that for large tan/3 there is a sizeable 
minimum value for the R-parity breaking decay branching ratio. 



9 Note that the MSSM LEP limits on neutralino masses do not hold in broken R-parity models. 



For a recent discussion see ref. [24]. 

10 These points correspond precisely to those we have eliminated from Fig. 2 as mentioned 
above. 



6 



4 Discussion 



We have shown that in a class of supersymmetric models with R-parity breaking 
the stop can have new decay modes to third generation fermions such as t\ — > b+r. 
For stop masses accessible at LEP2 and smaller than the chargino masses the 
branching ratio for the R-parity violating decay mode t± — > b+r can be comparable 
or even larger than the R-parity conserving mode t\ — > c + Xi- Note also that 
quite generally the R-parity breaking decay branching ratio could be increased 
through a larger value of h^s, a smaller value of 5, a larger value of the Higgsino 
content of the neutralino, or any combination of these. In fact, if R-parity is 
violated the lightest stop t± can be lighter than the lightest neutralino x± leading 
to a ~ 100 % branching ratio into the t± — > b + r decay mode. Finally, one should 
not forget that even in the case when ti — ► c+ Xi is the main stop decay mode, the 
resulting signature will be different from that expected in the MSSM, because 
the lightest neutralino will decay. From eq. ([!]), ©, © the expected modes are: 
Xi — > J + v, Sv, Z~ + /+ + u, r ± + l T + v, q + q + v, q + q' + v. While the first two 
decays give the same signature as in the MSSM, which is a c-jet and missing 
momentum, the others give rise to events with additional leptons or quarks. In 
order to identify the corresponding signals in detail, one would have to take into 
account all possible neutralino decay channels, and evaluate the corresponding 
detection efficiencies and standard model backgrounds in a dedicated Monte 
Carlo simulation. 

Acknowledgements 

This work was supported by DGICYT and Fonds zur Forderung der 
wissenschaftlichen Forschung, under grant numbers PB92-0084 and P10843- 
PHY, respectively. In addition, we have been supported by Accion Integrada 
Hispano-Austriaca (Biiro fur technisch-wissenschaftliche Zusammenarbeit des 
CAD", project no. 1), as well as by Doctoral Fellowships from DGICYT (M. A. G. 
J.) and CNPq-Brazil (M. B. M.). 



7 



References 



[1] H. P. Nilles, Phys. Rep. 110 (1984) 1; H. Haber and G. Kane, Phys. Rep. 
Ill (1985) 75 

[2] For a recent review see J. W. F. Valle, in Physics Beyond the Standard 
Model, lectures given at the VIII Jorge Andre Swieca Summer School 
(Rio de Janeiro, February 1995) and at V Taller Latinoamericano de 
Fenomenologia de las Interacciones Fundamentales (Puebla, Mexico, 
October 1995); |hep-ph/9603307 



[3] S. Dimopoulos, L.J. Hall, Phys. Lett. 207B (1988) 210; E. Ma, D. Ng, Phys. 
Rev. D41 (1990) 1005; V. Barger, G. F. Giudice, T. Han, Phys. Rev. D40 
(1989) 2987; T. Banks, Y. Grossman, E. Nardi, Y. Nir, Phys. Rev. D52 (1995) 
5319; M. Nowakowski, A. Pilaftsis, |hep-ph/9508"27T| ; G. Bhattacharyya, D. 



Choudhury, K. Sridhar, Phys. Lett. B349 (1995) 118; B. de Carlos, P. L. 
White, OUTP-96-01P 

[4] L. Hall, M. Suzuki, Nucl. Phys. B231 (1984) 419 

[5] C. S. Aulakh, R.N. Mohapatra, Phys. Lett. B119 (1982) 136; G. G. Ross, J. 
W. F. Valle, Phys. Lett. 151B (1985) 375; J. Ellis, G. Gelmini, C. Jarlskog, G. 
G. Ross, J. W. F. Valle, Phys. Lett. 150B (1985) 142; A. Santamaria, J.W.F. 
Valle, Phys. Lett. 195B (1987) 423; Phys. Rev. D39 (1989) 1780; Phys. Rev. 
Lett. 60 (1988) 397 

[6] A Masiero, J. W. F. Valle, Phys. Lett. B251 (1990) 273; J. C. Romao, C. A. 
Santos, and J. W. F. Valle, Phys. Lett. B288 (1992) 311 

[7] G. Giudice, A. Masiero, M. Pietroni, A. Riotto, Nucl. Phys. B396 (1993) 243; 
M. Shiraishi, I. Umemura, K. Yamamoto, Phys. Lett. B313 (1993) 89; see 
also I. Umemura, K. Yamamoto, Nucl. Phys. B423 (1994) 405 

[8] P. Nogueira, J. C. Romao, J. W. F. Valle, Phys. Lett. B251 (1990) 142; R. 
Barbieri, L. Hall, Phys. Lett. B238 (1990) 86; M. C. Gonzalez-Garcia, J W F 
Valle, Nucl. Phys. B355 (1991) 330; J. Romao, J. Rosiek and J. W. F. Valle, 
Phys. Lett. B351 (1995) 497; J. C. Romao, N. Rius, J. W. F. Valle, Nucl. Phys. 
B363 (1991) 369; J. C. Romao and J. W. F. Valle Phys. Lett. B272 (1991) 
436; Nucl. Phys. B381 (1992) 87 



8 



[9] J. W. F. Valle, Phys. Lett. B196 (1987) 157 

[10] V. Berezinskii and J.W.F. Valle, Phys. Lett. B318 (1993) 360; A. D. Dolgov, 
S. Pastor, and J.W.F. Valle frstro-ph/9506TJTT 



[11] J.-F. Grivaz, Rapporteur Talk, International Europhysics Conference on 
High Energy Physics, Brussels, 1995; ALEPH collaboration, CERN-PPE/96- 
10, submitted to Phys. Lett. B (1996) ;. H. Nowak and A. Sopczak, L3 Note 
1887, Jan. 1996 ; S. Asai and S. Komamiya, OPAL Physics Note PN-205, 
Feb. 1996 

[12] DO collaboration, Phys. Rev. Lett. 75 (1995) 618; and references therein, 
CDF collaboration, Phys. Rev. Lett. 69 (1992) 3439; UA2 collaboration, 
Phys. Lett. B235 (1990) 363 

[13] F. de Campos, M. A. Garcia- Jare no, A. Joshipura, J. Rosiek, J. W. F. Valle, 
Nucl. Phys. B451 (1995) 3-15; H. Hempfling, MPI-PhT/95-59 

[14] J. E. Kim, Phys. Rep. 150 (1987) 1; D. Dearborn, et al, Phys. Rev. Lett. 56 
(1986) 26; M. Fukugita et al, Phys. Rev. Lett. 48 (1982) 1522; Phys. Rev. 
D26 (1982) 1841; J. Ellis, K. Olive, Nucl. Phys. B223 (1983) 252 

[15] M. C. Gonzalez-Garaa, J. C. Romao, J. W. F. Valle, Nucl. Phys. B391 (1993) 
100 

[16] See, e.g. Fig. 25, ref. [|] 

[17] Proc. of the Workshop on Physics at LEP2, CERN 96-01, Vol. I, p. 463, eds. 
G. Altarelli, T. Sjostrand and F. Zwirner 

[18] A. Bartl, W. Majerotto and W. Porod, Z. Phys. C 64, 499 (1994), erratum 
Z. Phys. C 68, 518 (1995); A. Bartl, H. Eberl, S. Kraml, M. Majerotto and W. 
Porod, UWThPh-1996-18, HEPHY-PUB 642/96, [hep-ph/96034TO| , to appear 
in Z Phys. C 

[19] K. I. Hikasa and M. Kobayashi, Phys. Rev. D 36 (1987) 724 

[20] J. Ellis and D. V. Nanopoulos, Phys. Lett. B 110 (1982) 44; M. J. Duncan, 
Nucl. Phys. B221 (1983) 285; J. F. Donoghue, H. P. Nilles and D. Wyler, Phys. 
Lett. B 128 (1983) 55; M. Dugan, B. Grinstein and L. Hall, Nucl. Phys. B255 
(1985) 413; S. Bertolini, F. Borzumati, A. Masiero and G Ridolfi, Nucl. Phys. 
B353 (1991) 591 



9 



[21] L. J. Hall, V. A. Kostelecky and S. Raby, Nucl. Phys. B267 (1986) 415; 
F. Borzumati, A. Masiero, Phys. Rev. Lett. 57 (1986) 961; T. Kosmas, G. 
Leontaris, J. Vergados Phys. Lett. B219 (1989) 457; A. Gabbiani, A. Masiero, 
Nucl. Phys. B322 (1989) 235; R. Barbieri, L. Hall, A. Strumia, Nucl. Phys. 
B445 (1995) 219 

[22] Particle Data Group, Phys. Rev. D 45 (1994) 1173 

[23] A. Bartl, H. Fraas, W. Majerotto and N. Oshimo, Phys. Rev. D 40 (1989) 
1594 

[24] J.C. Romao, F. de Campos, M.A. Garcia-Jareno, M.B. Magro, J.W.F. Valle, 
|hep-ph/9604244 

[25] WA. Bardeen, M. Carena, S Pokorski and C.E.M. Wagner, Phys. Lett. B320 
(1994) 110 



10 



(a) 



(b) 



Figure 1: Feynman diagrams leading to the standard R-parity conserving decay 
(a) and the new R-parity violating stop decay (b). 




Figure 2: Barring unlikely small f\ values as discussed in the text, the curve 
displays the maximum values for the branching ratio for i\ — > b + r in per cent 
as a function of the lightest neutralino mass m^o. We have taken m~ tl = 80 GeV, 
Cl = 0.5, 5 = 0.1,v R = 100 GeV, v L = 0.1 GeV and h u33 = 0.03 and varied the 
remaining parameters as explained in the text 
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Figure 3: Scatter-plot for values of branching ratio for t\ — > b + r in per cent as a 
function of tan (3 for 55 < m xi < 60 GeV. All parameters are chosen as in figure 
2. 
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